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Abstract. We propose the reaHzation of linear crystals of cold ions which contain 
different atomic species for investigating quantum phase transitions and frustration 
effects in spin system beyond the commonly discussed case of s = i. Mutual spin- 
spin interactions between ions can be tailored via the Zeeman effect by applying 
oscillating magnetic fields with strong gradients. Further, collective vibrational modes 
in the mixed ion crystal can be used to enhance and to vary the strength of spin- 
spin interactions and even to switch those forces from a ferro- to an antifcrromagnctic 
character. We consider the behavior of the effective spin-spin couplings in an ion crystal 
of spin-i ions doped with high magnetic moment ions with spin S = 3. We analyze 
the ground state phase diagram and find regions with different spin orders including 
ferrimagnetic states. In the most simple non-trivial example we deal with a linear 
{Ca+, Mn+, Ca+} crystal with spins of {^, 3, ^}. To show the feasibility with current 
state-of-the-art experiments, we discuss how quantum phases might be detected using 
a collective Stern-Gerlach effect of the ion crystal and high resolution spectroscopy. 
Here, the state-dependent laser- induced fluorescence of the indicator spin-^ ion, of 
species *°Ca+, reveals also the spin state of the simulator spin-3 ions, ^°Mn+ as this 
does not possess suitable levels for optical excitation and detection. 
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1. Introduction 

Current ion trapping technology has led to rapid progress toward the realization of 
elementary quantum processors [1,2]. The ability to control the motional and internal 
states of the trapped ions with high accuracy allows for the experimental implementation 
of several textbook models such as quantum simulations of a Dirac equation with the 
Zitterbewegung and the Klein paradox [3,4]. On the other hand, the internal states of 
laser-cooled and trapped ions represent effective spins, which can be made to interact 
with each other for performing magnetic quantum phase simulations. These interactions 
may be realized by applying magnetic field gradients [5-10] or by laser light fields [11- 
15]. In both cases, spin-| ion crystals allow for a detailed investigation of complex 
quantum phase transitions and magnetic frustrated effects. Preliminary experimental 
steps in ion crystals have been realized with interacting spins of two ^^Mg+ [16] and 
three ^^^Yb+ ions [15, 17] in a linear configuration. Because of the long-range spin-spin 
interactions, the larger collections of trapped ions are expected to lead to intriguing, 
and so far unobserved phenomena, such as the formation of super-solids [18] or exotic 
quantum phases [19], where the specific advantages of the ion crystal are: almost perfect 
state preparation and readout with single site addressability, long coherence times, and 
a full tunability of the spin-spin interactions, even for long ranges beyond next-neighbor 
couplings. 

Going beyond spin-^ systems and trapping different ion species with spin S > \ 
will allow the study of novel aspects of quantum magnetism in a mixed spin chains 
[20,21]. Such impurity doped systems might model effects which are of interest in 
solid state physics [22,23]. Our proposal is inspired by the outstanding progress in 
quantum logic spectroscopy [24,25], where a single clock ion and a single readout ion 
are simultaneously confined and coupled through the mutual Coulomb repulsion, such 
that one can transfer the clock ion electronic state to the readout ion for high fidelity 
quantum state detection [26]. A different, new type of quantum logic readout technique 
enables us to propose quantum simulation in mixed ion crystals. For the case of neutral 
interacting atoms, the high magnetic moment of 6/xb of chromium has led to a wealth 
of novel effects [27, 28], made possible by the tuning of its spin interactions. ^°Mn, with 
an atomic number which is +1 higher as compared to Cr, shows a similar electronic 
structure and magnetic moment when singly ionized to Mn"*" for being trapped in the 
ion crystal. 

In this paper, we propose an efficient method for the creation of effective spin-spin 
interactions in ion crystals of spin-| ions doped with different ion species with spin S" = 3. 
An oscillating magnetic field gradient [8] can be used to implement coupling between 
the spin states of both ion species and the collective motional states of the impurity 
doped crystal, (see fig.l). The advantage to use an oscillating magnetic field instead 
of a laser field is due to avoiding the technical difficulties such as sideband cooling 
of the many vibrational modes and the necessity to use additional lasers to provide 
the spin-spin couplings. We show that by proper choice of the frequencies and the 
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Figure 1. Sketch of the proposed experimental sequence which a) initiaUzation in the 
static magnetic gradient field, followed by a b) RF tt/2 pulse on the Mn+ ions and 
the c) creation of the spin-spin interactions by using an oscillating magnetic gradient 
field, and finally concluded by the d) collective spin readout on the Ca+ ion using the 
position-dependent Stcrn-Gerlach effect. This leads to a bight or dark state of the Ca+ 
ion which is imaged on a CCD camera as state-dependent fluorescence. 

direction of the magnetic field gradient, the anisotropic Heisenberg model can be realized 
with tunable spin-spin couplings. We investigate the particular case of a field gradient 
applied along the trapping axis such that the spin-spin interactions are described by the 
transverse Ising model with single-ion anisotropy. We consider the ground state phase 
diagram for a small system consisting of two spin-| ions and one spin-3 ion placed 
at the center, which is realizable with the current ion-trap technology. Due to the 
complex competition between the spin-spin couplings and the single-ion anisotropy we 
distinguish four regions with different spin orders [15]. We find that for sufficiently strong 
antiferromagnetic nearest-neighbor coupling the spin order is ferrimagnetic wherein the 
two different spins are arranged in opposite directions. We show that the ferrimagnetic 
order can be frustrated due to competing next-nearest-neighbor coupling and the single- 
ion anisotropy which give rise to a highly entangled ground state. 

It has been shown that the atomic gases in optical lattice may be used to realize 
various condensed-matter models with high spin symmetry [21,30]. The ability to trap 
ions with large spins at a fixed position [31] and to tune the range and strength of the 
interactions makes the impurity doped ion crystal an analogue quantum simulator for 
quantum magnetism and frustration effects in a mixed spin system [32, 33]. 
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The paper is organized as follows: In sec. 2 we describe the theoretical background 
for the implementation of the effective spin-spin interactions in a [S, s) = (3, |) mixed 
spin system by using an oscillating magnetic field gradient. In sec. 3 we provide analysis 
of the ground state phase diagram for the transverse Ising Hamiltonian with single-ion 
anisotropy describing a spin system with two spin-| ions and one spin-3 ion placed at 
the center. The method for state preparation and readout of the spin states based on 
frequency addressing of the auxiliary ion is considered in sec. 4. Finally, in sec. 5 we 
give a summary of the results and discuss further and even more complex possibilities 
of quantum simulation with mixed ion crystals. 

2. Theoretical Model 

We consider a harmonically confined impurity doped ion crystal with N — K spin-| ions 
with mass m and K spin-3 ions with mass M. For instance, this is the case of ^^Ca"*" ion 
crystal doped with ^"^Mn"*" ions, which have ^5*3 electronic ground state. If the radial trap 
frequencies are much larger than the axial trap frequency {uJx,y ^ ^z)-, the ions arrange 
in a linear configuration along the axial z axis and occupy equilibrium positions [34]. 
The axial trap potential is independent of the mass, so that the equilibrium position of 
the ions are independent of the composition of the ion crystal. A static magnetic field 
Bq along the trap axis defines the quantization axis. The spin-| sublevels W) and \\) are 
Zeeman split by the applied magnetic field with a resonance frequency ujq = {gjfiB/h)Bo. 
Here gj denotes the Lande factor and /ig is the Bohr magneton. The spin sublevels of 
spin-3 ions possess seven Zeeman states which we index as |m) with magnetic quantum 
number m = —3, —2, . . . , +3 and resonance frequency cDq- In the case of ion crystal 
consisting of "^"^Ca"^ ions with electronic ground state ^6*1/2 doped with ^^Mn"*" ions the 
Lande gf-factor is (^j ~ 2 such that the resonance frequencies uq and ujq of both ion 
species are equal. 

2.1. Magnetic field gradient along the z-direction 

We assume that the impurity doped ion crystal interact collectively with oscillating 
magnetic gradient field with frequency u applied along the z direction (for simplicity 
we omit the constant magnetic offset) 

B = e^zEz cos ut. (1) 

The Hamiltonian for ions interacting with the magnetic field is H = Hq + Hj. Here 

hjjj ^ 

Ho = ^ ^ + ^o^Sl + Un,zai/an,z, (2) 

j=l k=l n=l 

is the interaction- free Hamiltonian, with being the Pauli matrix for the spin-| ion 
and S^. is the spin operator for the k^'^ spin-3 ion with Sz\m) = m\m). a^^ and an,z are, 
respectively, the creation and annihilation operators of collective phonons along the z 
axis with frequency UJn,z- The displacement zj of the j^^ ion from its equilibrium position 




Figure 2. left). The nearest-neighbor J12 = J23 = the next-nearest-neighbor J13 
spin-spin couplings and the single-ion anisotropy yl as a function of u for ion crystal 
consisting of two spin-i ions and one spin-3 ion placed at the center, (see Eq. (5)). The 
couplings J, Ji3 and A are normalized to e = {lS.zdzUiQY / {2uJz)-, which quantifies the 
change of the spin resonance frequency wq due to the shift of the equilibrium position 
of the ion with trap frequency lOz by an amount equal to the spread of the ground 
state wavefunction, Az = ^Jh/2rau}z- right). Quantum phases differ in the regions (I) 
to (IV). 



can be expressed in terms of these set of operators as Zj = "^^=1 ^j,n^^n{o,n^z + ^n,z)- 
Here Az„(a) = ^yh/2aun,z with a = m, M is the spread of the ground state wave 
function and (n = 1, 2, . . . , A^) are the normal mode eigenvectors in the z direction, 
(see Appendix A) [34,35]. The interaction between the magnetic dipole moment of the 
ion species and the magnetic gradient is described by Hj = —fl.B. The ^-component of 
the magnetic-dipole moment for spin-i ion is fi^ = {'y/2)az and, respectively, for spin-3 
fiz = 'jSz with 7 = fiBdL- We may transform the Hamiltonian in the interaction picture 
with respect to Hq to obtain 



N /N-K K \ 

n=l \ j=l k=l J 

N /N-K K \ 

E E i^^I + E ^InSl (at ,,e-(-— + a.,.e^(---)*) . (3) 

n=l \ j=l k=l J 

The function „ = bj ^AznB^'y /2h is the Rabi frequency of the j*'' ion, which quantifies 
the coupling to the n*^ vibrational mode. Hence, the oscillating magnetic gradient field 
mediates a coupling between the internal states of the ions and the external (motional) 
states of the ion crystal. Indeed, the two terms in the Hamiltonian (3) describe a time- 
varying spin-dependent displacement with frequencies (w -I- ujn,z) and (w — ujn^z)- If the 
frequency u is not resonant to any vibrational mode and the condition \uJn,z — ^ 
is satisfied for any n then we can perform time-averaging of the rapidly oscillating terms 
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in (3) [36]. Hence, we arrive at the following time-averaged effective Hamiltonian 

N-K K 

= l k.k' = l 

j>j' k>k' 

N K 
j,k=l k=l 

Therefore, the off-resonant oscillating magnetic gradient creates effective spin-spin 
interaction between the identical [14, 15] and different ion species in the crystal. The 
first two terms in (4) quantify the spin-spin coupling between the spin-| ions and the 
spin-3 ions. The third term in (4) describes the spin-spin coupling between the different 
ion species. Surprisingly, the adiabatic elimination of the vibrational modes for a ion 
crystal with s > | ions gives rise to single-ion anisotropy term Aj which quantifies 
the non-linear Zeeman shift of the spin-3 magnetic sublevels. We note that the non- 
linear Zeeman shift for spin-| ions is equal for the both magnetic sublevels, thereby it 
is proportional to the unit matrix. The couplings in (4) are given by 



n=l "'^ n=l "'^ 

.(3,2) _ B^-f Oj,n.Ofc,n .z _ B ^1 ST^ [Ok,n) /.s 



'n,z 



The main advantage of using an oscillating magnetic field gradient instead of constant 
is that we may engineer a variety of interactions between the ions. Fig. 2 shows the 
spin-spin couplings and the single- ion anisotropy (5) for a chain of two spin-| ions and 
one spin-3 ion place at the center versus the frequency u. In contrast with the constant 
magnetic field gradient applied along the trapping axis z wherein the spin-couplings can 
be only ferromagnetic, now as u is varied the magnitude and the sign of the couplings 
are changed which allows the creation of ferromagnetic, antiferromagnetic or frustrated 
interaction between the ions. The ground state of Hamiltonian (4) highly depends from 
the sign of the single- ion anisotropy terms A^. Indeed, for sufficiently large positive 
single- ion anisotropy {Al ^ 0), the spin-3 ions have magnetic quantum number m = 
for the ground state, while in the opposite limit {\Al\ ^ 0) the spin-3 ground state 
projection is m = ±3. 



2.2. Magnetic field along xyz-direction 

Consider the magnetic gradient applied along the xyz direction 

B{x,y,z) = e^zB^cosut + e^xB^f{t) - eyyByf{t), (6) 

with f{t) = (coscjftt + cos Urt). Such a field can be created in a micro structured planar 
ion trap, which contains a central wire loop [9]. The oscillating field in x-y direction 
provides additional coupling = 6j_^Ag„i?q7/2/i (g = x,y) between the internal 
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Figure 3. The ground state phase diagram calculated by an exact diagonalization 
of Hamiltonian (8). a) The ferromagnetic population Pf function of the 

normalized transverse magnetic field B^/e and the frequency uj. In region I) the 
spin-spin couplings are ferromagnetic J, J13 < and the single-ion anisotropy is 
^ < 0. b) The population difference Pf^i — + Pf,3 — Pa,3 as a function of the 
normalized transverse magnetic field B^/e and the frequency w. In region II) the 
spin-spin couplings are antiferromagnetic J, J13 > and the single-ion anisotropy 
is A > 0. By increasing w as i?° — the system undergoes transition \ijjf\i) — )■ 
\iljf,2) — > |'0/,3) ~^ IV'a,3)- c) The ferrimagnctic population as a function of the 
normalized transverse magnetic field B^/e and the frequency uj. In region III) the spin- 
spin couplings are, respectively, J > 0, J13 < and the single-ion anisotropy is A < 0. 
d) The antiferromagnetic population Pa_o as a function of the normalized transverse 
magnetic field B^/e and the frequency lo. In region IV) the spin-spin couplings are 
J < 0, Ji3 > 0, respectively, and the single-ion anisotropy is A > 0. 

and motional degree of freedom of the ion crystal. We assume that the frequencies 
Ub — ujq = 6 and —5 are tuned to the blue- and red- sideband transition 

with detuning ±6. Then the x-y gradients induce a spin-dependent displacement with 
frequencies {6 + C0n,q) and {6 — uJn,q) similar to Eq. (3), where ujn,q are the vibrational 
frequencies, respectively, in the x and y directions, (see Appendix A). After applying 
an optical rotating-wave approximation (neglecting the terms cjq + tUfe.r) and assuming 
that \ujn,q — 5| ^ is fulfilled for any vibrational mode in the x-y-direction, the 
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time-averaged effective Hamiltonian is given by 

N-K K 



q=x,y,z jj'^i k,k'=i 
j>j' 

TV K 



j,k=l k=l 

The spin-spin couplings in x-y direction are identical in form to (5) by replacing u ^ 5 
and z ^ x,y. Hence, the magnetic field (6) creates an anisotropic Heisenberg (XYZ) 
interaction between the effective spins. 



3. Transverse Ising Model 

The quantum transverse Ising Hamiltonian is given by 

N-K X K 

Hti = - hBl{ E Y + E ^^'}- (8) 

3=1 k=l 

The last term in (8) can be simulated by driving transitions between the ion spin states 
employing radio frequency field Bq = CxBq cos ut. Assuming the resonance condition is 
fulfilled, i.e., u = uq 'we obtain the effective transverse field B^ = 'yBo/2h. 

The simplest non-trivial case is to consider ion chain with two spin-i ions and 
one spin-3 ion placed at the center. Such ordering of ions is consistent with the 
natural behavior, as observed in Ref. [31]. When applying the oscillating gradient 
field, the resulting spin-spin couplings are shown in fig. 2 as a function of the drive 
frequency u. We may distinguish four different regions wherein the spin-spin interactions 
are ferromagnetic, antiferromagnetic, or frustrated. The presence of the single-ion 
anisotropy A in (8) changes substantially the ground state phase diagram compared 
to the case of spin-| ion chain. In contrast to the spin-| string, the ground state of 
the mixed {S, s) = (3, |) spin system can be frustrated due to the complex competition 
between the spin-spin couplings J and J13 and the single- ion anisotropy A. 

In region I), see fig. 2, all interactions are ferromagnetic ( J, J13 < 0) and the 
resulting ground state of the Hamiltonian (8) as i?^ — )■ is a coherent superposition 
of two ferromagnetic states |tt) |3) and \ID |— 3). Because the single-ion anisotropy is 
negative {A < 0), the ground state energy is minimized for spin-3 state with magnetic 
quantum number m = ±3. The ferromagnetic population -P/m,3 = -Pfts + -^U-s 
a function of the effective magnetic field and the frequency u is shown in fig. 3a. 
In the region II) the spin interactions are antiferromagnetic ( J, J13 > 0). This is 
the case of ferrimagnetism in which spins of two types interact by nearest-neighbor 
antiferromagnetic coupling, (see fig. 3b). However, the ferrimagnetic interaction 
is frustrated due to the competing next-nearest-neighbor antiferromagnetic coupling 
Ji3 > which disturbs the ferrimagnetic order and tends to align the two spins- 1 in 
a antiferromagnetic state. Additionally, the ferrimagnetic interaction is also frustrated 
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due to the strong positive single-ion anisotropy A > which attempts to project the 
spin-3 state |m) in a quantum number m = 0. In the beginning of region IT) A > 
is high but the ferrimagnetic configuration is still energetically favorable such that the 
ground state of Hamiltonian (8) as -B^ — ?■ is a superposition of two ferrimagnetic 
states \i>f,i) = (Itt) + |l))/v^! wherein the spins S = 3 and s = ^ are aligned 
anti-parallel with each other. By increasing u the single-ion anisotropy A decreases 
which allows the quantum number m to increase and the resulting ground state is 
\i^2,f) = (Itt) |-2) + |ii) |2))/V2. For A/J = 2/3 which occurs at w ^ 1.24a;^ the ground 
state for 5° = is an entangled superposition of four ferrimagnetic states |tt) 
Itt) |1), Itt) l~2), and |tt) |2). By decreasing A the ferrimagnetic state \ip2,f) undergoes 
a transition to \^3j) = (|tt) l~3) + |tt) |3))/"\/2- At the transition point u = lAuz and 
A/ J = 0.4 the resulting ground state for 5^ = is an entangled superposition of 
four ferrimagnetic states, |tt) |-2), |tt) |2), |tt) |-3), and |tt) |3). For sufficiently high 
positive Ji3 the ferrimagnetic arrangement is not favorable and the two spin-i ions 
are arranged in an antiferromagnetic state. Hence, the ferrimagnetic configuration is 
broken and the ground state is iV^a.s) = (Itt) + ltt))(|3) + |— 3))/2. At the transition 
point u ~ 1.65Uz, the ground state of the Hamiltonian (8) is an entangled superposition 
of six states: two ferrimagnetic states |tt) l~3) and |tt) |3), and four antiferromagnetic 
states Itt) |3), |tt) |-3), |tt) |3), and |tt) |-3). In region HI) the nearest-neighbor 
spin coupling is antiferromagnetic (J > 0) such that the ferrimagnetic configuration 
is favorable. In contrast with the region II), no frustration exist because the next- 
nearest-neighbor coupling is ferromagnetic J13 < and A < 0. In the entire region the 
ground state is a ferrimagnetic IV^/,3), which minimizes all interactions. In fig. 3c the 
ferrimagnetic population Pf^^ = P|f_3 + P^s versus the effective magnetic field and the 
frequency u is shown . In region IV) the single-ion anisotropy is positive {A > 0) which 
causes the spin-3 being projected into quantum number m = 0. Additionally, the next- 
nearest-neighbor coupling is antiferromagnetic J13 > and the resulting ground state 
as 5° — )■ is antiferromagnetic \ipa,o) = (Itt) + Itt)) |0) /V^- The antiferromagnetic 
population Pa^ = P-^io + P^to is shown in fig. 3d. 

3.1. Towards larger and more complex mixed crystal 

By increasing the number of ions the spin-spin interactions become more complex and 
consequently the new ground states exhibit a variety of spin orders. In alternating 
mixed-spin chain we can find regions in which the nearest-neighbor interactions are 
antiferromagnetic with ferrimagnetic ground states, (see fig. 4). Such antiferromagnetic 
mixed-spin chains have been proposed as a model for describing certain molecular- 
based magnets of experimental interest [37]. Moreover, one of the most studied topics 
in alternating mixed-spin chains is the appearance of a plateaux in the magnetization 
curve [33, 38-40]. It has been shown that the formation of these plateaux depend highly 
on the competing interactions and single-ion anisotropy. Thus, the proposed ion trap 
based simulator of quantum magnetism in mixed-spin system will allow for detailed 
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Figure 4. a). The spin-spin couplings and single-ion anisotropy in alternating 
mixed-spin chain {^,3, ^,3, ^} consisting of three spin-i ions and two spin-3 ions. 
All couplings are normalized to e. The frequency of the oscillating magnetic field 
is set to cj = l.Sujz- The nearest-neighbor couplings are antiferromagnetic. Similar 
antiferromagnetic behavior is observed in larger alternating odd number ion crystals, 
b) Because of the antiferromagnetic nearest-neighbor interactions, the different spins 
are aligned anti-parallel with each other and the resulting ground state as B^/e — ?> 
is ferrimagnetic |Vj/,2> = (Itt) |-2,-2) + |U> |2,2))/v^ or |^/,3) = (Itt) |-3,-3) + 
14,4,) |3,3))/-\/2. The figure shows the population difference — as a function of 
w and B^/e. 



studies of such phenomena. 

4. Preparation and Spin Readout 

In order to prepare the initial state of the Hamihonian 

< = -^5°{f + ^2+f} (9) 

and to readout the final state after the adiabatic ramping of the spin couplings we 
can apply a static magnetic field gradient B = Bq + bze^ along the trapping axis 
z which shifts the resonance frequencies of transitions between Zeeman states, (see 
fig. 5a and the general scheme in fig. 1). Due to the spatial variation of the magnetic 
field, the spin states exhibit site-specific resonance frequency, u'^ = uq + 'jbz/h. Due 
to the applied gradient field, the equilibrium position of each of the ions is shifted 
depending of the spin states of the other ions in the linear crystal. Indeed, the total 
force which acts on the magnetic dipole moment of spin-| ion along the trapping axis 
^ is Fi = fP + fP, where fP = hd,uj{S,) and fP = {h/2)d,oo{a,) are forces 
associated with the spins ^5*2 and {h/2)az in a magnetic field B and () indicates 
the expectation value [5]. The spin-dependent force shifts the equilibrium position 
by an amount = {Fz^^^ + Fz^"^) / {mul). Consequently, the spin resonance frequency is 
changed from uji to u'^, with uj'i= uji + 'ybdz/h, (see fig. 5b). The frequency separation 
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Figure 5. a) In the presence of static magnetic field gradient the Zeeman shift of the 
spin states is different for each ion. b) The ion equihbrium position is shifted due to 
the spin-dependent force acting on the magnetic dipole moment. The direction and 
the magnitude of the displacement depends on the spin states of the other ions. Thus, 
the spin resonance frequency will change from ui to • By scanning the frequency of 
the applied laser light the spin transition comes into resonance and consequently the 
ion scatter light. 



SuJzeeman = — ^i\ between the two spin resonances of spin-| ion is given by 

S^Zeeman = ^ ^ \dz\- (10) 

Assuming a high but reahstic magnetic gradient of 6 = 20 T/m [10] and an axial trap 
frequency of = 27r x 100 kHz, the shift from the equihbrium position of a spin-i ion 
corresponding to the state \m = 8)2 |t)3 is approximately ~ 50 nm which gives rise 
to frequency separation (10) of about Suzeeman ~ 174 kHz. Note that much smaller 
frequency shifts are easily resolved on a narrow Raman or quadrupole transition, such 
that in a final step the observation of laser-induced fluorescence allows one to determine 
the spin states. The experimental sequence is started by preparing the two spin-i ions in 
a state \ID l^'s), where l^^s) is a still unknown state of the spin-3 ion. Additionally, we 
may trap an auxiliary ion initially prepared by optical pumping in the state which 
can be used to readout the spin states of the remaining three ions. In the second step we 
switch on the static magnetic gradient, which creates the spin-dependent force. Next, 
we expose the auxiliary ion to laser light with frequency uj[, which drives the transition 
li)a ^ lt)a oiily if the state of spin-3 ion is |m = —3). If the auxiliary ion does not 
scatter light we have to discard the measurement and restart. If we observe fluorescence 
then the state of the spin-3 ion is measured in |— 3) and results, after performing a 7r/2 
rotation along the y axis, is the desired ground state |\E'o) of Hamiltonian (9) 

1^0) = e-'f'^ie-'5^2e-'5'^' |-3) . (11) 
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In the same way the probabihty observable Ps^^s3,m. for the state jsiSa) |m) can 
be determined. By scanning the laser frequency, the auxiliary spin-| ion comes into 
resonance and consequently scatter light. The auxiliary resonance frequency is defined 
depending of the spin states of the other three ions. The readout method cannot be 
applied directly for a general collective spin state with more than one spin-3 ion. For 
example, in a mixed spin system with two spin-^ ions and two spin-3 ions, the states 
Itt) \m, —m) and Iff) |— m, m) create equal amounts of displacement to the equilibrium 
position of the auxiliary ion, such that both states are indistinguishable. For these 
more general cases, the state detection problem can be solved by (i) applying a spatially 
varying magnetic field gradient with d'^B/dz^ ^ which allows one to distinguish the 
spin states due to the position dependence of the force, or (ii) a separating the entire ion 
crystal into smaller sub-crystal and transporting them into magnetic gradient detection 
zones [41, 42] which is a promising solution especially for segmented micro traps [9]. The 
very precise measurement of the position wave function of an ion was an experimental 
key element in experiments which could demonstrate the quantum random walk [43], 
was also used in other ion trap experiments with single ions and relies on the light ion 
interaction in either Raman [44] or bichromatic light fields [3,43]. These methods can 
determine the wave packet position to an accuracy of a few percent of the ground state 
wavepacket extension A^;^, but for an application here one would need to lock the laser 
light phases for position determination with respect to the initialization RF pulses on 
the ions, which has not been experimentally realized so far. 

5. Conclusion and Outlook 

We have proposed a method for the creation and manipulation of the spin-spin 
interactions in spin-^ ion crystal doped with high magnetic moment ions with spin 5 = 3. 
It is shown that by tuning the frequency and direction of the oscillating magnetic field 
gradients various fundamental models in quantum magnetism of mixed spin systems can 
be realized. Because of the competing long-range spin-spin couplings the spin orders 
are extremely numerous even for spin system consisting of a small number of ions. We 
have proposed a technique for spin preparation and readout based on the frequency 
addressing of an indicator spin-| ion in the presence of spatially varying magnetic field. 

In future we will investigate the decoherence properties of such states, and how the 
quantum simulation is affected by the typical noise sources in an ion trap experiment. 
As the different quantum states are affected very differently by the expected dominating 
sources of noise, the ambient magnetic field fluctuations, we expect that states in 
decoherence-free subspaces, as for example \il)f^i) = (|tt) + oi' I'V'a.s) = 

(Iti) + l4-t))(|3) + l~3))/2 with zero effective magnetic moment are preferred. A second 
future research theme is the combination of spin-dependent forces by magnetic gradients 
and with an advanced control of multiple segments in a micro trap [9, 45] which form the 
axial potential for the confinement of the ion crystal in the z-direction. As this potential 
may be shaped in non-trivial, non-harmonic way and even with multiple potential wells 
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[46,47] the interspatial distances in the ion crystal become free parameters instead of 
being given by the Coulomb repulsion in an overall harmonic potential. Here, it was 
shown that 2x N cluster states can be generated [48] as a resource for one way quantum 
computing. 
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Appendix A. Calculation of the eigenfrequencies and eigenvectors of the 
impurity doped ion crystal 

In order to determine the spin-spin interactions we need to calculate the eigenfrequencies 
and eigenvectors of the impurity doped ion crystal. The potential energy of collection 
oi N — K ions with mass m and K ions with mass M is 

, so N-K , so N-K 

i=i i=i 



9 N 



(A.l) 



where for simplicity of the notation we assume that the trapping confinements in the 
radial x-y direction are equal. Here cJz(m) ~ m and Uz{M) ~ M 

are, respectively, the axial trap frequency for ion with mass m and M. The radial 
oscillation frequency uj^.{a) ~ with a = m, M is reduced by the axial trap frequency 
according to 



V 2a;r°(a)'/ 

in lowest-order approximation [49]. For small displacement the motional degrees of 
freedom in x, and z directions are decoupled. The Hessian matrices Aij and Bij 
which describe the small oscillation of the ions around their equilibrium position in the 
axial z, and respectively, x-y directions are given by 
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Here we have introduced the mass ratio fi = M/m and anisotropy parameter a = 
u z{rn) / uj^{m) and Uj being the dimensionless equihbrium position of the j^^ ion. In 
the above expressions the index jm denotes the position of the j^^ ion with mass M. 
The eigenfrequencies Uz^n = ^z{jn)\/Ki and eigenvectors {n = 1, 2, . . . , A^) in the z 
direction are given by the diagonahzation of the following matrix 

j 1^ hi): 

= 3 = Jm), 



such that Y.f=i ^ij^ln 



Aj/fJ', ( 

Ajl^, (z or J =JM, «7^j)> 
Aij I II, {iu ^ Jm), 

= XnMn- same way the eigenfrequencies Ug^r. 



(A.5) 



w°(m) 



and eigenvectors (g = x, y) in the radial direction are obtained by the diagonahzation 
of the following matrix 

Bij, (i = j ^ 3m), 
Bij/fi, {i =j = Jm): 
B^, {i^3^3m): (A.6) 
Bij/ y/J^, {i or j = Jm, i 7^ j), 
Bij/fi, {in 7^ Jm), 



Bij 



with Ell Bijbl^ = -inhln 
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